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Abstract 

We give a concise proof of nonlinear stability for a large class of solutions to the Einstein equations 
with a positive cosmological constant and compact spatial topology, where the spatial metric is Einstein 
with either positive or negative Einstein constant. The proof uses the CMC Einstein flow and stability 
follows by an energy argument. We prove in addition that the development of non-CMC initial data 
close to the background contains a CMC hypersurface, which in turn implies that stability holds for 
arbitrary perturbations. Furthermore, we construct a one-parameter family of initial data such that 
above a critical parameter value the corresponding development is future and past incomplete. 


1 Introduction 

Understanding the long time behavior of the Einstein flow and the global geometry of the resulting space- 
times has been a major field of interest in General Relativity for the past 30 years. A particularly success¬ 
ful area concerns the nonlinear stability problem for explicit solutions to Einstein’s field equations. The 
first general results are due to Friedrich IIFr86l for deSitter space-time and Christodoulou-Klainerman 
IChK1931 for the Minkowski space-time. Since then several results of a similar nature for different back¬ 
grounds have been established. In this paper we focus on the Einstein-flow with a positive cosmological 
constant. 


1.1 The positive cosmological constant 


The late time asymptotics of homogeneous cosmological models in the presence of a cosmological con¬ 
stant have been first analyzed by Wald lWa83 1. Following Friedrich’s work on the stability problem 
HFr861 [Fr86-1 1 for the 3+1-dimensional case, Anderson generalized the stabilty result to asymptotically 
de Sitter space of arbitrary even dimension | An05 1. Eater, Ringstrom was able to find conditions on the 
initial data, such that the global evolution problem can be localized to a coordinate neighborhood HRiQ8H . 
Using this local result he showed stability and future completeness for large classes of initial data on 
arbitrary spatial topologies. This implies that the spatial topology itself cannot be deduced from the long 
time behavior of the Einstein flow in the presence of a positive cosmological constant. A result of this 
nature had also been established by Friedrich in HFr86-11 . Ringstrom’s results hold for the Einstein-scalar 
field system and he later generalized them to the Einstein-Vlasov-scalar field system in HRi 131 . Similar 
results have been obtained by Svedberg for the Einstein-Maxwell system BSvl 1 1 . 

In the case of the Einstein-Euler system questions of long-time existence are complicated by the likely 
appearance of shocks. However, in the presence of a positive cosmological constant it has been shown 
by Rodnianski and Speck that the accelerated expansion is sufficiently strong to avoid shock formation in 
the non-vacuum setting for the irrotational Einstein-Euler system [ RoSpl3] and by Speck for the general 
Einstein-Euler system [Spl21. 
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1.2 The CMC-Einstein Flow on compact manifolds 

In the study of nonlinear stability of expanding solutions to the vacuum Einstein-flow with vanishing 
cosmological constant CMC (constant mean curvature) foliations have been proven to be very beneficial. 
The study of the CMC Einstein flow was initiated by the work of Andersson, Moncrief and Tromba on 
the global existence of CMC foliations of vacuum solutions of the Einstein equations in 2+1 dimensions 
!AnMoTr971 . Fischer and Moncrief I FiMoO 11 IFiMo02l studied the Einstein flow in CMC gauge for the 
higher dimensional case, which eventually led to the proof of stability for FLRW(Friedmann-Lemaitre- 
Robertson-Walker) type solutions in 3+1 dimensions by Andersson and Moncrief l AnMo()4 | and finally 
to the stability of a large class of spatial Einstein geometries of negative scalar curvature in arbitrary 
dimensions by the same authors in 2011 HAnMol 11 . The proof in l AnMol 1 1 is based on a carefully 
adjusted energy argument which shows the asymptotic convergence of the perturbed solution to the spatial 
Einstein metric. 

The motivation for the present paper is the study of the CMC Einstein flow with a positive cosmological 
constant, which has so far only been considered for the 2+1 dimensional case by Andersson, Moncrief 
and Tromba in IAnMoTr971 . 

1.3 Main Results 

The aim of this paper is to prove nonlinear stability for the Einstein flow in CMCSH gauge ( constant 
mean curvature spatial harmonic gauge ) with a positive cosmological constant for arbitrary spatial di¬ 
mensions n > 2 and different background geometries. 

We begin by constructing homogeneous model solutions for the CMC-Einstein flow with positive cos¬ 
mological constant which have spatial Einstein geometries with positive and negative Einstein constants. 
In particular, we construct a one-parameter family of initial data, such that there is a critical parameter 
value above which the corresponding future and past development recollapses, while for smaller values 
the developement is future and past complete. This family illustrates the variety of scenarios which might 
occur despite the presence of a cosmological constant. Our main result is the proof of nonlinear stabil¬ 
ity for a large class of spatial Einstein-geometries with positive as well as negative Einstein constants in 
arbitrary dimensions n > 2. In addition, we show that the future development of arbitrary perturbations, 
not necessarily those which solve the CMC constraints, is geodesically complete towards the future and 
admits a CMC-foliation. The main theorem in the case of a negative Einstein constant is the following. 


Theorem 1.1. Let M be a smooth compact n-dimensional manifold (n > 2) without boundary and 7 be 
an Einstein metric satisfying Ric( 7 ) = — (n — 1 ) 7 . Then for s > n/2 + 2, s' > n/2 + s and e > 0 there 
exists a <5(e) > 0 s.t.for initial data (go, fco) satisfying 


llso 


fc 0 + V2"i 


< 5 

H s'-l 


( 1 . 1 ) 


its maximal globally hyperbolic development under the Einstein flow with positive cosmological constant 
A = can be foliated by CMC-hypersurfaces M t , t £ [arcsinh(l), 00 ) such that the induced 

metrics gt satisfy 

||sinh - 2 (t)< 7 t - 7 || ffs < e. ( 1 . 2 ) 

In particular, all corresponding homogeneous solutions are orbitally stable and the future developments 
of small perturbations are future geodesically complete. 


In the case of positive Einstein constant, the main theorem is the following. 

Theorem 1.2. Let M be a smooth compact n-dimensional manifold (n > 2) without boundary and 7 be 
an Einstein metric satisfying Ric( 7 ) = (n — 1)7 which does not admit Killing vector fields and such that 
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—2 (n — 1) is not an eigenvalue of the Laplacian. Then for s > nil + 2 , s' > n/2 + s and e > 0 there 
exists a (5(e) > 0 s.t. for initial data (go, ko ) satisfying 

llso - l\\ H °' + PolliP'-i < s d- 3 ) 

its maximal globally hyperbolic development under the Einstein flow with positive cosmological constant 
A = n ^ n f 1 l can p e globally foliated by CMC-hypersurfaces M t , f £ K such that the induced metrics g t 
satisfy 

||cosh -2 (£)g t - 'y\\ HS < e. (1.4) 

In particular, all corresponding homogeneous solutions are orbitally stable and the future- and past 
developments of small perturbations are future- and past geodesically complete, respectively. 

The idea of the proof is based on an energy argument, which makes use of the elliptic hyperbolic struc¬ 
ture of the Einstein flow in CMCSH gauge and which is inspired by the ideas of Andersson and Moncrief 
in l AnMol 1 11. The presence of the positive cosmological constant yields a specific asymptotic hierarchy 
of the terms appearing in the CMCSH equations, which needs to be taken into account by choosing an 
appropriate rescaling of the evolving geometry. The eventual energy estimate does not contain a decay 
inducing negative term on the right had side as exploited by a correction mechanism in BAnMolll . Here, 
another crucial observation allows to obtain a sufficient energy estimate in the small data setting. The 
idea is as follows. The energy is essentially the sum of a geometric Sobolev norm of the difference be¬ 
tween the rescaled spatial metric and the Einstein background geometry and that of a geometric Sobolev 
norm of the trace free part E of the rescaled second fundamental form. The Sobolev norms are defined 
w.r.t. a Laplace-type operator on tensors corresponding to the background geometry and the perturbed 
spatial metric. Straightforward energy estimates for this norm contain a negative term on the right had 
side, which is a multiple of the Sobolev norm of E. This term on its own cannot be exploited to gain ad¬ 
ditional decay for the energy - also a correction mechanism fails in this case. However, a careful analysis 
of all additional perturbation terms on the right hand side of the energy estimate allows to use this term 
to absorb all perturbation terms with insufficient decay properties under the assumption that we are in a 
small data scenario. We obtain a strong energy estimate with an exponentially decaying coefficient under 
the smallness condition. Global existence and stability then follow by a bootstrap argument and, taking 
the rescaling into account, the asymptotics of the solution imply the desired completeness. In addition, 
we demonstrate that the nonlinear stability result also applies to non-CMC perturbations. In particular, 
we prove that the maximal globally hyperbolic development of any initial data close to the background 
contains a CMC surface, to which we can apply our analysis and obtain future completeness and a corre¬ 
sponding CMC foliation of these spacetimes. 

The main theorem which we prove contains results of previous papers mentioned above. In IIRi08l 
Ringstrom considers the case of positive curvature in dimensions n > 3 and that of negative curva¬ 
ture in dimension 3. However, for these solutions the existence of CMC foliations has not been proven 
to our knowledge. In addition, we do not have to exclude the 2 dimensional case. In the case of positive 
curvature our assumptions imply n > 3, since we have to exclude the case of the sphere. The reason why 
our approach fails for the sphere is probably due to the expected non-existence for CMC foliations for 
certain initial data close to the sphere. 


1.4 Overview of the paper 

In section[2]we introduce our notations and recall the equations for the Einstein flow in CMCSH gauge 
with positive cosmological constant. In subsection 2.4 we construct homogeneous solutions for the cases 
of positive and negative spatial Einstein manifolds. In subsection 2.5 we prove the existence of CMC 
surfaces in the MGHD (maximal globally hyperbolic development) of non-CMC initial data and also the 
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generality of the spatial harmonicity condition. Section [3] contains a discussion of the reduced Hamilto¬ 
nian for the Einstein flow with positive cosmological constant analogous to BFiMo021 . In particular we 
obtain similar monotonicity results as IFiMo02l for general initial data and we note them for the sake of 
completeness. They are not essential for the remainder of the paper. In section[4]we construct the family 
of initial data containing both data with expanding and recollapsing asymptotics. Section [5] contains the 
proof of the main theorem. In subsection |5.2| we address the problem of local existence for the data we 
consider. The proof of the main theorem is divided into the analysis of the elliptic system in subsection 
|5.3| the main energy estimate in subsection |5.4| an improved decay estimate in subsection |5.5| and the 
bootstrap argument which establishes global existence in subsection [576] 
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1.5 Notations and Conventions 


We collect all relevant further notations used in this paper in the following. We define the scalar product 
of two tensors u, v w.r.t. 7 by 

(u,v) = u^vu Y kr y jl - (1-5) 

In addition, we define the mixed L 2 -scalar product by 


( u > v )L 2 (g, 7 ) 



( 1 . 6 ) 


The corresponding L 2 -scalar product where also the volume form is taken w.r.t. 7 is denoted by (. , .) 7 . 
We denote the standard Sobolev spaces on M by H s , where the norm is defined w.r.t. a fixed metric 
on M. We do not distiguish the notation for Sobolev spaces of different types of tensors. We denote 
the iT s -norm by INI The ball in the H s topology of radius e > 0 around a tensor h is denoted by 
B s e {h). The covariant derivative of a certain metric g is denoted by V [g\. If it is clear from the context, 
the reference to the metric is suppressed. The curvature tensor is defined with the sign convention such 
that RijkiX k = Vj V, X[ — V, VjV/. The Laplacians are defined with the sign convention such that all 
eigenvalues are nonpositive. 


2 The CMC-Einstein flow 

2.1 Background solutions 

Throughout the paper, we put the cosmological constant to A = n ( n ~ 1 ) where n is the spatial dimension. 
Then Einsteins equations are equivalent to RiCg = ng. Assume that (M, 7 ) is a compact Riemannian 
Einstein manifold with negative scalar curvature R( 7 ) = —n[n — 1). Then, the metric 

7 = —dt 2 + sinh 2 (f )7 ( 2 . 1 ) 

solves Einsteins equation on (0, 00 ) x M. This solution is future geodesically complete and the mean 
curvature of the slices M t = {t} x M is r(f) = — n which is strictly monotonically increasing 

on (0, 00 ) and tends from —00 as t —> 0 to —n as t 00 . Introducing r as a new time variable, 7 is 
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a solution of the CMCSH flow we will introduce in the next section. Similarly, if 7 is Einstein but of 
positive scalar curvature R( 7 ) = n{n — 1), 

7 = — dt 2 + cosh 2 (f) 7 ( 2 . 2 ) 

solves the Einstein equation on R x M which is future and past geodesically complete. If (M, 7 ) is 
the sphere, we recover the de-Sitter metric and therefore these models are called generalized de-Sitter 
spaces. The mean curvature of M t is r(f) = —n which is strictly monotonically decreasing. We 
have lim t _ J . 00 r(f) = — n and lim t _ > _ 00 r(f) = n. The metric 7 cannot be regarded as a solution of the 
CMCSH flow, but as a solution of the reversed CMCSH flow if we put the time variable to —r. 

Remark 2.1. If 7 is a Ricci-flat metric, 7 = —dt 2 +e 2t 7 also solves the Einstein equation. Unfortunately, 
the mean curvature of any slice M t is constantly —n and therefore, we cannot handle this solution with 
our gauge conditions. 

2.2 ADM Einstein equations 

We consider a space-time of the form ffi. x M, where M is a smooth compact /(-dimensional manifold 
without boundary. For the Lorentzian metric we choose the ADM-Ansatz 

(rl+1) <? = — N 2 dt ® dt + gijiydx 1 + X l dt) (g) ( dx J + X j dt), (2.3) 

where g = ( N , X, g) denote the lapse function, shift vector field and the spatial metric, respectively. The 
Einstein equations in CMCSH gauge. 


tl 'gk = T = t 

w'dt/ - f %) = v k = 0, 


(2.4) 


where T*-, I’d denote the Christoffel symbols w.r.t <j and 7 , respectively, with positive cosmological 
constant A = ’’dn-i) reac | 


R(g)~\n 2 +r 2 


n — 1 


= n(n — 1) 


V l E ij = 0 


d/fZ / 'y — 2 . V (E j j -)- T/ng i:/ j A C x gij 

dt^ij = N(Rij + t Ejj — 2EjfeE j + (r 2 /n — n)gij ) 

+ 23xE ij - gij -Ej j — ViVyjA 

n n 


E -| n 

y n 


AN = -1 + A 

AX i + R l m X m - C x V i = + r(2/n - 1)V 4 A 


- (2iVE m " - {C x g)™){r mn - r mn ) 


(2.5) 


where the second fundamental form k has been decomposed into 

k = E + -g, (2.6) 

n 

where E denotes the tracefree part. Note it is assumed that N > 0. We have used the following standard 
notations. R(g) denotes the Ricci scalar curvature of g, C x denotes the Lie derivative w.r.t. the shift, Rij 
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denotes the Ricci tensor of the metric g. The Laplacian A is understood to be defined w.r.t. g. 

Finally, we remark that in the case of an reversed CMC-gauge, t = —r, which we use for spatial Einstein 
metrics of positive curvature, one has the lapse equation in the form 


AN = 1 + N 




(2.7) 


The equation for the trace free part of the second fundamental form in this case reads 


dt^ij = N(Rij + rT,ij - 2EtfEj* + (r 2 /n - n)gij) 
1 2Nt 

+ Cx^ij H —gij S ij — ViVjN 

n n 

and the other equations remain the same. 


( 2 . 8 ) 


2.3 The Einstein operator 

The fundamental property of the CMCSH-Einstein flow lies in its elliptic-hyperbolic structure given by 
the decomposition of the spatial Ricci tensor into the Einstein-operator and perturbation terms as given 
in the following lemma. 

Lemma 2.2 ( ll AnMolTi Lemma 6.2]). Let 7 be an Einstein metric, then we have the expansion 


R( 7 ) 1 

Rij 5ij g^ — 2 ^ 9 , 1(9 l)ij + Jij: 

(2.9) 

where 

O 

£g,-yh = —A g ^h — 2R 1 h 

is the Einstein operator. The Laplacian has the local formula 

(2.10) 

A g^hij = — V[7] TO (5 m? VgV[7]„/iij), 

Tg 

(2.11) 

O 

the curvature action is given by R 7 (h)ij = R t k ji ( 7 ) h ki for some 2-tensor h, Sij = 
J is an error term which satisfies the estimate 

\ (Vj Vj + Vj Vi ) and 

\\J\\H-i<C\\g-'Y\\ 2 H .- 

(2.12) 


2.4 Homogeneous solutions of the CMC-Einstein flow 

We recover the background solutions discussed before in the present gauge by assuming homogeneity. In 
the homogeneous setting, meaning a vanishing trace free part of the second fundamental form, a spatially 
constant lapse function and vanishing shift vector we obtain the following solutions. 

Standard CMC-gauge 

For the standard CMC-gauge, t = t, we deduce from the lapse equation that 

Ti 

N = - 2 - ^ (2.1.3) 

t z — 

and since N > 0, r 2 > n 2 . Then the evolving physical metric is given by 

9{j) = 9{t 0 ) 1 \—~2 j (2.14) 

T- — 
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where the metrics have the property that the Ricci tensor is given by R^ = — IL jf g l: j. We recover the 
metric written in CMC time. In particular, the scalar curvature of the Einstein manifold is negative. 


Reversed CMC-gauge 

Suppose now, we have the reversed CMC-gauge, i.e. t = —r. Then 

Ti 

N = ~2 - 2 ( 2 . 1 5 ) 

n z — r z 

and therefore, t 2 < n 2 . The solution for curve of the physical metric is the same as above, where 
Rij = gij, i.e. the scalar curvature is positive. We recover \2.2) . 

2.5 Universality of the gauge conditions 

An important issue arising in the context of the CMCSH gauge concerns the generality of perturbations 
which can be evolved by the CMCSH Einstein flow. Considering CMC initial data induced by a back¬ 
ground solution which admits a CMC foliation one would prefer to consider general perturbations of this 
initial data, i.e. solutions to the general constraint equations without the CMC condition. But such non- 
CMC initial data cannot be evolved by the CMCSH Einstein flow. 

This problem can be overcome by the following construction. Assuming non-CMC initial data close to 
the initial data induced by the background solution, which is CMC, general theory assures the existence 
of a maximal globally hyperbolic development of this initial data. Under the smallness condition it is 
possible to show that this development contains a CMC surface. Starting the evolution from this CMC 
surface one can analyze the geometry of the corresponding future development and eventually treat all 
perturbations by considering the equations in CMCSH gauge. 

A similar question arises for the spatial harmonic gauge, which also is shown to apply to general initial 
data. We discuss these aspects of the gauges in the following. 

2.5.1 Spatial Harmonicity 

Let Ai be the set of smooth Riemannian metrics on M. Fix a metric 7 £ A4 and let 77, be the set of 
metrics g £ A 4 such that id : (M, g) — > (M , 7) is a harmonic map. In other words, 

n = {g G M I V g k = g ij (T[g% - r[ 7 ]* ) = 0} , (2.16) 

where T [ 5 ], T [7] are the Christoffel symbols of g, 7 , respectively. Our aim in this section is to prove that 
under certain conditions on the background metric 7, R is a smooth submanifold of AA and a local slice 
of the action of the diffeomorphism group through 7 . By the first variation of the Christoffel symbols 
(see e.g. !Be08 . Theorem 1.174]), the differential of the map <I> : < 71 —>• is given by 

d$(hy = V j hj - ^Vhigh, (2.17) 

where h is a symmetric 2-tensor and V is the covariant derivative w.r.t. g. 

Lemma 2.3. Let (M, 7) be an Einstein manifold such that — 2 /n ■ R( 7) is not an eigenvalue of the 
Laplacian A 7 and 7 does not admit Killing vector fields. Then the operator 


P:X l A 7 X i + i?J[ 7 ]X J ' 


( 2 . 18 ) 
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is an isomorphism which presents the decomposition 

X(M) = {grad/ | / G C°°(M)} © {A G C°°{TM ) | divA = 0} . (2.19) 

Proof. We suppress the dependance on 7 in the following notation. By a standard argument using com¬ 
mutators of covariant derivatives, we have 

AV7 + Ky j f = V*A/ + 2iRV7 = V‘Af + — ■ V*/ (2.20) 

j j n 

which shows that because of the eigenvalue assumption, P maps the first factor bijectively onto itself. By 
self-adjointness of P, the second factor is also preserved. We define maps L and L* by 

L-.X^^ViXj+VjXi), L* :h^r -X 3 hj. (2.21) 

Note that L* is the adjoint map of L with respect to the L 2 -scalar product induced by 7 . Now for any 
vector field A' with divA = 0, we have 

(. L*LXy = V j V j X i + V j V%- = AX 1 + X j X l Xj - V l X 3 X j = AX 1 + R)X 3 = (PA) 1 . 

( 2 . 22 ) 

Thus, PX = 0 implies LX = 0. But the kernel of L contains precisely the Killing vector fields, and 
hence. A' = 0. Therefore, P is injective and by self-adjointness, P is also suijective. □ 

Lemma 2.4. Let (M, 7 ) be an Einstein manifold such that —2/n ■ R( 7 ) is not an eigenvalue of the 
Laplacian A 7 and 7 does not admit Killing vector fields. Then, d4> 7 : C°° (S 2 M) —> C°°(TM) is 
surjective. Moreover, we have the splitting 

C°°(S 2 M) = ker(d$ 7 ) © im(L). (2.23) 

Here, S 2 M denotes the bundle of symmetric 2-tensors. 

Proof. To prove the first assertion, we consider a vector field A and compute 

d $ 7 o L(Xy = ^(V J VjA* + V J V i Aj) - |rv,\' j = AA‘ + R) X 3 = {PX)\ (2.24) 

Due to Lemma [23] P is an isomorphism, so <i<1 > 7 is surjective even when restricted to Lie derivatives. 
This calculation also shows that ker(V/< 1 > 7 ) n im(L) = 0. To prove that the direct sum spans all of 
C°°(S 2 M), we first note that 

ker(d<£> 7 ) © im(d$*) = C°°(S 2 M) (2.25) 

so it suffices to prove that 

im(d<I> 7 ) C ker(d4> 7 ) © im(L). (2.26) 

If h G im(d'f , 7 ), there is a vector field A such that 

h = d$;(A) = -LX + idiv(A) • 7 . (2.27) 

We may use an arbitrary vector Y field to rewrite this expression as 

h = —L(X + Y) + LY + ^div(A) • 7 (2.28) 
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We are done with the proof if we can find Y such that 


LY + ^div(X) • 7 £ ker(d<h 7 ). 


(2.29) 


Thus, we have to solve 

0 = d<f> 7 (LF + *div(AT) • 7 ) = PY + *d$ 7 (div(X) ■ 7 ) (2.30) 

where we used \2.2A) . This can be done due to Lemma [23] □ 

For our purposes, it is more convenient to work on neighbourhoods with Sobolev regularity. We therefore 
use H s -norms with s > ” + I for the following theorem. We remark that the above lemmas also hold, 
if we descend to H s -regularity. Let M s be the space of // '-metrics on M and let 7 ~L S be the set of all 
g £ M s satisfying the condition in ( |2. 1 6[ >. 

Theorem 2.5. Let (M, 7 ) be an Einstein manifold such that — 2/n ■ R( 7 ) is not an eigenvalue of the 
Laplacian A 7 and 7 does not admit Killing vector fields. Then in a small H s -neighbourhood IA C M s 
of 7 , Tl s is a smooth submanifold of M s with tangent space 

T^n s = j/r € H s (S 2 M) I = Y j h) - ^VTr h = oj . (2.31) 

Moreover, for any g £lA there exists an isometric metric g £ P s which is H s -close to 7 , i.e. there exists 
ip £ Lf s (Diff(Af)) such that g = p*g. 

Proof The first assertion follows from the first assertion of Lemma [2~4] The second assertion follows 
from the implicit function theorem for Banach manifolds applied to the map 

V : U s x 7T s (DifF(M)) -> M s (2.32) 

given by ip) = p*g. Since there are no Killing fields, d v I ; ( 7 .i c j) is injective and its image is 

im(dtf (7|id) ) = © {Cxi I X £ H S (TM)} = her(d$ 7 ) © im(L) (2.33) 

which equals H S (S 2 M) = 1{A4 S by ( | 2.23 [ >. Therefore, is a diffeomorphism from a /-/ '-neighbourhood 
of ( 7 , id) in Li s x FP(DifT(M)) to a ^“-neighbourhood of 7 in M s . □ 

Remark 2.6. The assertions of Theorem |2. 5 1 hold for any Riemannian metric 7 where the operator P is 
an isomorphism. 

2.5.2 Constant mean curvature hypersurfaces 

Let us now consider the CMC-gauge. Let M be a compact manifold, JcMan open and bounded interval 
and M k ’ a the set of C A: '' t -Lorentz metrics on I x M, such that the induced metrics on the hypersurfaces 
Mt = {t} x M are all Riemannian. A Banach manifold structure on this set is induced by the norm 

m = H^lUo + + \\9t\\c^ , (2.34) 

where we identify g according to the foliation by submanifolds M t with the triple (/V, A', g t ) of the Lapse 
function, the shift vector and the induced metrics g t = g\ivr t ■ The norms of the right hand side are taken 
with respect to the Riemannian metric dt 2 + 7 on / x M. Let C k ’ a {M 1 I ) the set of C k,a -functions 
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/ : M —> I endowed with the natural Banach manifold structure. Each such function defines naturally 
an embedding : M —» I x M by if{x) = (f(x), x ). We define a map 

H : M k ’ a x C k+1,a (M, I) D V -» C k ~ 1,a (M) (2.35) 

which associates to each pair (g, f) the mean curvature along the embedding if : M —> I x M induced 
by the metric g. Here, V is the open subset of pais (g, /) such that (■ if)* 9 is a Riemannian metric. 

This is a smooth map between Banach manifolds. To see this, it suffices to consider a local ex¬ 
pression of this map: Using local coordinates on M, we see that Trf( x \ x )im(if) is spanned by the 
vectors dif(di) = dif • dt + di. Let F(t,x) = t — f(x) and the matrix (gf)ij be defined by ( g/)ij = 
g(dif(di),dif(dj)). Then the mean curvature is 

H SJ = v = (2-36) 

where v is the timelike unit normal and (<jfY J is the inverse of ( gf)ij . This expression contains second 
derivatives of the function / and first derivatives of g. We use an implicit function theorem applied to the 
map H to prove the following lemma. 

Lemma 2.7. Let (M, 7 ) be a compact Einstein manifold with scalar cun’ature R( 7 ) = —n{n — 1) and 
let I be an arbitrary open and bounded interval around 0. Let £>\ and consider the metric 

7 = —dt 2 + cosh 2 (t) 7 

whose initial data induced on the hypersurface {0} x M is ( 7 , 0 ). Then for any C e,a X C^~ 1,a - 
neighbourhood IA of { 7 , 0 ), there exists a neighbourhood V C of 7 such that any g £ V admits 

a hypersurface such that the pair ( g, k ) of the metric and the second fundamental form induced on this 
hypersurface is in U and ti g k = 0 . 

Proof. Consider the map II of above and note that H ( 7 ,0) = 0. We compute its differential at the tupel 
( 7 ,0) restricted to the second argument. By the variational formula of the mean curvature in ll BaBrCaOS . 
Proposition 2.2], 

dHz[fl(f),w) = — [A 7 w — (Ri&y(<9 t , d t ) + |fc| 2 )iu] = — [A 7 uj + nw\. (2.37) 

Because we excluded the case of the sphere, the operator A 7 + n : C f ' +1,a {M) —> C^ 1,a (M) is 
an isomorphism | Oh62 Theorem 1 and Theorem 2], Due to the implicit function theorem for Banach 
manifolds, we have neighbourhoods U' C C t,a {M) of 7 , V' C C t+} ' a (M, I) of 0 and a smooth function 
F : W —> V such that H(g,F(g)) = 0, i.e. F associates to each metric g a minimal Riemannian 
hypersurface given by the graph of the function F(g). Moreover, F(jj) is the only function in V' such 
that H(g, F(g)) = 0. The proof is finished by the remark that the map g —► (< 7 , k ) associating to g the 
metric and the second fundamental form of graph(F(^)) is continuous from V to C e,a x C' £_1,a . □ 

Lemma 2.8. Let {M, 7 ) be an Einstein manifold with scalar cun’ature R( 7 ) = — n(n — 1) and let I be 
an arbitrary open and bounded inten’al in (0, 00 ) around £0 = arcsinh(l). Let £ > 1 and consider the 
metric 

7 = —dt 2 + sinh 2 (£)7 

whose initial data induced on the hypersurface {arcsinh(l)} x M is ( 7 , —\/2 r )). Then for any C e,a x 
C ?- 1 "-neighbourhoodU of ( 7 , —s/Qr/), there exists a neighbourhood V C C e,a (f4) ofy such that any 
g £ V admits a hypersurface such that the pair (g,k) of the metric and the second fundamental form 
induced on this hypersurface is in IA and tr g k = —\/2n. 
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Proof. The proof is analogous as above. In this case, we consider the map H = H + \[2n. We have 
H(j, to) = 0 and we compute 

dHg tto ( 0,tu) = — [A 7 u> — (Ric g(d t ,d t ) + |fe| 7 )iu] = —[A 7 ru — nw\. (2.38) 

which is always an isomorphism from C l+1,a (M) to C e ~ 1,a (M). □ 

Theorem 2.9. Let 7 be an Einstein manifold of scalar cunmture R( 7 ) = n(n — 1) (resp. R( 7 ) = 
—n(n — 1)) and let s > n /2 + 1, s' > n/2 + s. Then for any H s x H s—1 -neighbourhood U 3 ( 7 , 0) 
(resp. U 3 ( 7 , — s/2^)) of CMC initial data sets, there exists a H s x H s -neighbourhood V 3 ( 7 , 0) 
(resp. U 3 ( 7 , —\f2y)) of general initial data sets such that any development of initial data in V admits 
a CMC-hypersurface such that the initial data induced on the hypersurface lies in U. 

Proof. Let ( g i: kf) be an initial data set converging to ( 7 , 0) in H s x H s -1 . By the proof of | Ri09 , 
Theorem 15.10], one obtains a sequence of solutions of Einsteins equations gt such that for each slice 
{1} x M, t £ I, the data (JVj, X^, ( gt ) t ) converges in H s to the corresponding data of the background 
solution 7 . Moreover, we have H s _1 -convergence of the time derivatives (ffN,. dfXf) by the choice of 
the gauge used in the proof of the above mentioned theorem. Moreover, we have H s _1 -convergence of 
d t gi to the corresponding quantity of 7 . 

By Sobolev embedding, we obtain convergence of ( Ni,Xi , (gf)t) in C s,a and convergence of their 
time-derivatives (d t Ni, dtXi, dtgi) in C s ~ 1,a on each slice {f} x M. Using the gauge condition and the 
Einstein equation, we also obtain convergence of higher time-derivatives of the above quantities so that g, 
converges to 7 in C s,a {I x M). For i large enough, the metrics g t admit hypersurfaces of constant mean 
curvature due to the lemmas above and the initial datas (g, . kg) induced on the hypersurfaces converge in 
C s ’ a x C S_1 ’ Q , hence also in H s x H s ~ x . This proves the theorem. □ 


3 Reduced Hamiltonian 

The behavior of the reduced Hamiltonian under the Einstein flow with vanishing cosmological constant 
has been studied by Fischer and Moncrief in lFiMo02l| . They show that the reduced Hamiltonian is 
monotonically decreasing in expanding direction of the flow and its infimum equals the topological sigma 
constant of the spatial manifold. In the following section we present some analog results for the Einstein 
flow with positive cosmological constant for the sake of completeness. Results of this type are clearly 
interesting due to the fact that they hold for arbitrary, not necessarily small data. They might therefore be 
relevant in a study of large initial data, which is in a reasonable sense far from the background solutions. 
In the present paper, this section is independent of the remainder and might therefore be skipped without 
affecting the understanding of the proof of nonlinear stability. 


3.1 Monotonicity and stationary points 

Theorem 3.1. Let g t be a solution of the CMC Einstein flow with positive cosmological constant A = 
which is expanding (i.e. r < 0). Then the reduced Hamiltonian 


( T 2 V 72 

(— -nj vol St , tg(-oo,- n) 


(3.1) 


is monotonically decreasing and stays constant if and only if g t is a family of (up to isometry) homothetic 
metrics which are Einstein metrics of negative scalar curvature. 
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Proof. The time-derivative of the reduced Hamiltonian is 

T (vol gt -C—~ n ) [ ■ (3- 2 ) 

\ n Jm J 

By integrating the equation on the Lapse function, we can replace the right hand side and we get 

/ 2 \n/2-l 

dtH~ d =[-~ n ) T m\ 2 P 9t > 0 (3.3) 

which proves the first part of the theorem. Suppose now that the Hamiltonian stays constant. Then E = 0, 
since N is positive. Since E = 0 holds on any time interval, where the reduced Hamiltonian is constant, 
we also have d t E = 0 on this interval. Since the operator A — (—— n) is invertible, the equation on 
the lapse function is uniquely solvable. Thus, N = (—— n) -1 . If we now insert E = 9 t E = 0 in the 
evolution equation for E, we see that the equation Rij = -' holds for all time. It follows 

from the evolution equation of the metric, that the g t are homothetic up to isometry. □ 

Theorem 3.2. Let g t be a solution of the reversed CMC Einstein flow with positive cosmological constant 
A = which is expanding (i.e. r < 0). Then the reduced Hamiltonian 

( t 2 \ n/2 

H+ed = - — J v ol St , t e (-n, -n) (3.4) 

is monotonically increasing. If H+ ed is constant along the flow, gt is a family of constant positive scalar 
curvature metrics which lie up to diffeomorphism in the same conformal class. Moreover, if—R(g t )/(n — 
1) is not an eigenvalue of A gt , g t is a family of homothetic metrics which are Einstein metrics of positive 
scalar curvature. 

Proof. The time derivative is 

d tH+d = d -rH+ed = ( n ^ (~ vol st + - ff) J N Pa^ ■ (3-5) 

By integrating we can replace the right hand side by 

dtH+ d =(n-~) ’ (-r) / N\T,\ 2 p gt > 0, (3.6) 

V n J JM 

which proves the monotonicity. As above, E = 9 t E t = 0 as long as the Hamiltonian stays constant. By 
the Hamiltonian constraint and the evolution equation of the metric, 

R{g) = (n- l)(n- —), d t g = -2N-g + C x g, (3.7) 

n n 

which proves the first assertion in the equality case. If —R(g t )/{n — 1) is not an eigenvalue of A gt , the 
equation on the lapse is uniquely solvable and so, N = (n — By the evolution equation for E, this 

immediately yields R^ = N~ lli= l-gij. As in the previous theorem, the metrics gt are homothetic up to 
isometry in this case. □ 


dtH rP d — — 


— n 


i/2— 
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3.2 Critical points on the reduced phase space 

Since these Hamiltonians are (up to the time-dependant scale factor) the same as in the case of vanishing 
cosmological constant IFiMo02l . one can also do the same analysis as in this paper. Suppose the manifold 
M is of negative Yamabe-type, i.e. the scalar curvature of any metric on M is negative somewhere. 
Then in any conformal class, there is a unique metric of constant scalar curvature —1. By the standard 
conformal method, we may therefore consider the reduced phase space 

v; ed = { (g, E) G M x T(S 2 M) \ R(g) = -1, tr s E = 0, (V^Ey = 0} . (3.8) 

If we drop the assumption of negative Yamabe-type, then the uniqueness statement about metrics of fixed 
scalar curvature in a conformal class fails. On the other hand, given a metric g close enough to an Einstein 
metric of constant scalar curvature 1 which is not the round sphere, then there exists a unique metric of 
scalar curvature 1 in the conformal class of g ltBoWaZi04l Theorem C]. Therefore, we may in this case 
introduce the phase space 

V+ d = {( ff ,E) e M x T(S 2 M) | R(g) = l,tr ff E = 0, (V 9 )%- = 0} . (3.9) 

One may now regard the reduced Hamiltonians H± d as functionals on the Phase spaces V^r ed and develop 
their variational theory as in HFiMo021 Section 4]. The critical points are the pairs (g , E) where E = 0 
and Ric g = ± —g, respectively. The Hessian of H~ ed at a critical point will be positive semidefinite if 
and only if the Einstein operator is positive semidefinite on TT-tensors (transverse and traceless 
tensors) whereas the Hessian of 11+, d at critical points will be negative semidefinite if and only if H gg is 
positive semidefinite on TT-tensors. In the case of negative Yamabe-type, the global infimum of H~ ed , 
will furthermore determine the topological sigma constant which is also known as the Yamabe invariant 
(c.f. [FiMoOZ Theorem 6 ]). 


4 One-parameter family of initial data with collapsing and expand¬ 
ing regimes 

This section is concerned with the construction of a one-parameter family of initial data such that for a 
parameter value strictly above a certain threshold, the future and past development recollapses while for 
the critical value and below the corresponding future and past development expands for all time. The 
initial data consists of a product of positive Einstein metrics with identical Einstein constants while both 
metrics are multiplied by a large respectively by a small constant - yielding non-equilibrium initial data, 
where the small factor recollapses. For parameter values close to 1 the Einstein metrics are initially of 
almost similar volume and expand both for infinite time. The initial data corresponding to the threshold 
value of the parameter yields a solution where one factor remains constant in time while the second metric 
expands for infinite time. We proceed with the explicit construction. 


We consider a product manifold M x N such that g,\f and g\; be Einstein metrics of positive scalar 
curvature on M and N, respectively, with dim M = dim N = m. Let n = 2m. The Einstein constants 
are chosen such that 

Ric 9N = (n- 1 )g N and Ric 9M = (n - 1 )g M . (4.1) 

Given s € ( 5 , 00 ), let 

§ 

9 m ( s ) = s • g M and g N (s) = — — • ffjv• (4.2) 

zs — 1 

We consider now a Lorentzian metric 

9 = -dt 2 + a(t) 2 g M {s) + b{t) 2 g N (s) (4.3) 
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on I x M x N, where / C Si is some interval, g is supposed to be a solution of the Einstein equations 


Ricg = ng 


(4.4) 


with initial conditions 

a(0) = 6(0) = 1 and a'(0) = 6'(0) = 0, (4.5) 

which are compatible with the constraints. Furthermore, we define the new variables, x = logo and 
y = log 6 . Einstein equations then imply the system of ODE’s 

x" = n — -(n — l)e~ 2x — — [(a/) 2 + x'y'] 

S i n ( 4 -6) 

y" = n - (2 - -)(n - l)e 2y - ~{y' 2 + x'y'} 


with initial data x(0) = y(0) = x'(0) = y'{ 0) = 0. The equation f?oo = n<?oo yields 


n . a" b" 

2 ( tt + y ) = " 


(4.7) 


and equivalently 


x" + y" + x' 2 + y' 2 = 2 . 


If s = 1, we recover the generalized de-Sitter metric since a(t) = b(t) = cosh(f) in this case. 


(4.8) 


Theorem 4.1. Let s £ (|,oo). Consider initial data (<7m(s)© <7at(s),0) on M x N where gM{s), gN(s) 
are as in ( |4.1| > and ( |4.2| >. Then, for n f} 2 < s or s < , the future and past development is geodesically 

incomplete. For s £ t] the future and past development is geodesically complete. Moreover, if 

s £ ( 1 , n ~i), we have a limit 

v n ’ n—2 


(7, = lim 


vol (M,a(t) 2 g M (s)) 
t->±oo vo\(N 1 b(t) 2 gN(s)) 


(4.9) 


Proof. Let us prove the first assertion. Without loss of generality, we restrict to the case s > ■ By 

Lemma 4.2 below, there exists a time To > 0 such that lim t ^r + y(t) = —oo which in turn implies that 
the scale factor of the metric gjy(s) satisfies lim t _>,T + 6 (f) = 0. By choice of the initial values, a and 6 are 
time-symmetric, hence lim t ^._x + 6 (f) = 0 as well. Clearly, these solutions are geodesically incomplete 
in the future and the past. 

Let us now prove the second assertion. The case s = 1 is the case of th e de -Sitter space, so there is 
nothing to prove. We may restrict to the case 1 < s < Then by Lemma 4.3 the functions x, y hence 
the scale factors a, 6 exist for all f > 0. By time-symmetry, they exist for alTT 


£ JR. All these solutions 
are future- and past geodesically complete. Finally, the existence of ( |4.9| ) follows from ( |4.2()[ >. 

□ 


4.1 Evolution of non-equilibrium initial data 

We consider first the case of non-equilibrium initial data, which we define by 

n — 1 

s > - -z. (4.10) 

n— 2 

If s < R —- the roles of pm{s) and pn(s) interchange, so we may restrict to the first case. We prove the 
following lemma. 
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Lemma 4.2. For non-equilibrium initial data initial data, the solution (x(t),y(t)) of the system ( |4.6| > 
with initial data x(0) = y( 0) = a/(0) = j/'(0) = 0 does not exist for all time. More precisely, there exists 
a time T + > 0 such that lim t _ > 7 ’ + y(f) = — oo. 

Proof The condition on s implies 

y"(0) =n- (2- -)(n- 1) < 0 (4.11) 

5 

and 

x"(0) =n- -(n- 1) > 2. (4.12) 

s 

We show that x and y are strictly monotonically increasing (decreasing, respectively) on the interval 
of existence. 

Since a/'(0 ) > 0, we have x' > 0 for small t > 0. Let to > 0 denote the first time, such that 
x'(to) = 0. Then ( |4.6| ) implies 

x"{t Q ) =n~-(n- l)e- 2x(to) > 2 > 0, (4.13) 

s 


which in turn implies x'(t) < 0 for t £ (fi — e, t±). Therefore x' > 0 as long as it exists. Analogously 
one can show that y is strictly monotonically decreasing on the interval of existence. Next we show, that 
x exists at least as long as y" und y'. By ( |4.8[ > and the monotonicity of x and y we obtain 

- y" = (n - 1)[(2 - -)e~ 2y ^ e ~ 2x )] ™[(z') 2 - (l/) 2 ] 

s s 2 (4. 14 ) 

< 2(n — 1 ) — ^[(x') 2 — {y 1 ) 2 ] 

Using ( |4.8| >, x" can be eliminated from this inequality. Elementary manipulations yield 

Or ') 2 <4 +^_ 2 /'+^( 2 /) 2 . (4.15) 

which proves the claim. We show in the following that y is unbounded from below. Assume the contrary. 
Then by the strict monotonicity of y the existence of a limit lim^^ y(t) =: y{oo) follows, so does 
limt^oo y'(t) = 0. Then 

/•OO 

/ y"{t)dt = lim y'(t) - y'{ 0) = 0, (4.16) 

Jo *“*- 00 

so either lim^oo y"(t) = 0 or there is a sequence s* — > oo such that y”(si ) = 0 for all i £ R Using 
( |4.6| > would then imply x'(t) diverges for t —> oo, which however is a contradiction to ( |4. 1 5[ i. 

Finally, we show that y blows up in finite time. Addition of both equations in ( |4.6[ ) yields 

x" + y" = 2n - -(n - l)e" 2x - (2 - -)(n - l)e~ 2y - ^(x' + y') 2 . (4.17) 

s s 2 

By monotonicity of x and y and the unboundedness of y from below there is a constant C > 0 and a time 
1 1 > 0 such that the differential inequality 


+ y" <-C-^{x' + y') 2 


(4.18) 


holds for all t>t\. Here t\ > 0 is arbitrary. The corresponding ODE is solved by the tangens. Therefore 
x' + y' blows up in finite time (say t^ > tf) towards — oo. As x' > 0, y' blows up. In addition, we also 
have y(tf) = J* 2 y'(t)dt = — oo, which implies that y diverges. □ 
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4.2 Equilibrium initial data 

We consider now equilibrium initial data given by 

n — 1 

l<s< -(4.19) 

n — 2 

Lemma 4.3. For equilibrium initial data initial data, the solution (x(t),y(t)) of the system d4~6> with 
initial data tc(0) = y( 0) = x'(0) = y'( 0) = 0 exists for all t > 0. Moreover, if s < we have a limit 

C s = lim (x(t) - y(t)). (4.20) 

t—> OO 

Proof. We consider the case of 1 < s < Then 

y , 7 (0) = n — (2 — -)(n — 1) € (0,1) (4.21) 

S 

and 

*"(0) = n-^(n-l)e (1,2). (4.22) 

Both x and y are strictly monotonically increasing. We have x'(t) > 0 for small t. Let to be the first time, 
where x'{to) = 0. Then by ( |4.8| > we obtain 

x"{to) = n — -(n — i) e ^ 2 x ( t o) > o, (4.23) 

s 

as long as y'(to) exists. Thus x' < 0 on (to — e, to), which causes the contradiction. Analogously one 
shows that y is strictly monotonically increasing. From \A.\1\ and mo no tonicity we deduce 

2 - + y') 2 < x" + y" <2 n-^(x' + y') 2 . (4.24) 

The solution of the corresponding ODE is tanh. This implies 0 < x' + y' < C for all t > 0. Due to the 
positivity of x' and y 1 these statements hold for x' and y' individually. In particular, x and y exist for all 
times. In addition we have x'(t ) + y'{t) > C\ > 0 for all t > t\ and x(t) + y(t) > ^t for all t >t\ 
and C 2 > 0. Using what we have seen so far we obtain the following estimates. 

x" ^ y" = (n - 1)[(2 - -)e~ 2y -e~ 2x )] - ™l(x') 2 ( y ’) 2 } 

s s 2 

< (n - 1)(2 - l )e- 2 ^ x+v) - ™{x' - y')(x' + y') (4.25) 

5 2 

<(n^l)(2~-)e- 2C ^-C 3 (x'-y') 

s 

This differential inequality holds for t > t\. On the other hand, we also have 

x" — y" > —(n — l)-e~ 2C2t — C 4 (x' — y'). (4.26) 

s 

From these inequalities we deduce that x l — y' decays exponentially and converges to 0 as t —t 00 . The 
exponential decay implies the existence of the limit 

/»oo 

lim (x(t) - y(t)) = / {x’it) - y\t))dt. (4.27) 

t~> °° Jo 
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In the boundary case s = y"(0) = 0. Thus, y = 0 and the system reduces to the initial value 

problem 

x" = n — {n ~ 2)e~ 2x — ^(a/) 2 , x(0) = x'(0) = 0. (4.28) 

By similar arguments as above, one shows that x is strictly monotonically increasing. An immediate 
implication is x" < n, which implies that x grows at most quadratic ally. Therefore, it it exists for all 
time. □ 


4.3 Products of negative curvature 

Finally, we address the case of products of negative Einstein metrics. We consider an analogous con¬ 
struction as in the previous sections. Let qm and be two compact ?n-dimensional Einstein metrics 
with 

RiCjm =-{n- 1 )g M and Ric gjv = ~(n~l)g N , (4.29) 

where n = 2m. Let s £ (|, oo) and 

§ 

9 m{s) = s ■ g M 9 n{s) = - - - ■ gN (4.30) 

Zs — 1 

We consider a Lorentzian metric of the form 

9 = —dt 2 + a(t) 2 g M {s) + b(t) 2 g N (s ), (4.31) 

and demand Ric g = n ■ g. We have the conditions a( 0) = 6(0) = 1 and a'(0) = 6'(0) = y/2 compatible 
with the constraints. Defining the variables x = log(a) und y = log( 6 ) yields the system of ODE’s 

x" = n + -(n — l)e~ 2x — ^[{x') 2 + x'y'] 

s 2 (4.32) 

y" = n+ (2- -){n- l)e~ 2y - ^[( y') 2 +x'y'] 

with initial conditions x(0) = y( 0) = 0, x'(0) = y'{ 0) = \[2. Due to f?oo = n ■ g oo. we have 


n.a" b" 

2<^ + T )= " 


or equivalently 


x" + y" + (x 1 ) 2 + {y') 2 = 2 . 

In the case s = 1, we recover the background metric 

— dt 2 + sinh 2 (f — arcsinh(l))(^M © ffjv)- 


(4.33) 

(4.34) 

(4.35) 


For this system we obtain the following result. 


Theorem 4.4. Let s £ (|,oo). Consider initial data ( 3 m(s)® ffjv(s), ^\/ 2 ( 3 M(s)® 5 Ar(s))) on M x N 
where gutis). gN{s) are as in ( |4.29| ) and ( |4.30| >. Then, the future development is geodesically complete. 
Moreover, we have a limit 


r = y vol(M, a(t) 2 g M (s)) 

3 ti™ vol(N,b(t) 2 g N {s)) ‘ 


(4.36) 


Proof. For any s £ (|, oo), an analogoue of Lemma 4.3 
theorem follows as in the second part of Theorem 4.1 


can be proven by the same arguments. Then the 

□ 
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5 Nonlinear Stability 


We turn now to the main part of the paper, presenting the proof of the nonlinear stability results. Theorem 
1 1.1 1 and Theorem |1.2| The proof consists of four steps: a rescaling of the system, local stability results, 
elliptic estimates for lapse and shift and a uniform energy estimate for the evolving data g and £. The steps 
distinguish formally between the cases of positive and negative curvature of the background geometry. 
However, both cases can be handled more or less similarly, so that we present most of the arguments only 
for one case explicitly and in detail. 


5.1 CMC-Einstein flow in rescaled time 

Let us relabel the solution (g. E, B, X) of the elliptic-hyperbolic system by g. E, N, X. We now intro¬ 
duce new scale-invariant variables by 

9 = s(r)g, N = s(t)N, £ = ^t) 1 / 2 ^ X = s(t) 1 ' 2 X. (5.1) 


5.1.1 Rescaling of the CMC flow 

In the CMC case, we define the scale factor as s(r) = (-) 2 — 1. In this variables, the constraint equations 
read 

R(g) — |E| 2 = — (n — l)n, V Q E ab = 0. (5.2) 

Furthermore we define a new time-variable T by the equation r = — . This time coincides with 

the time of the solution ( |2.1| ). We now rewrite this system in the rescaled variables and the time variable 
T. The defining equations for lapse and shift are 

A g A = -l + A(|£| 2 +n) 

A g X z + R m X m = 2 VjN&i - cosh(T)(2 - n)X l N - (2iV£ mn - {C x g) mn ){r mn - f ! mn ). 

(5.3) 

Here we additionally used that VjE*- 7 = 0 and V 1 = g ,J (T* : ; — f Jj ) = 0 The evolution equations are 


„ cosh(T) . . 

drg%j = - 2 —— — ((1 - nN)gij - 


sinh(T) 
,cosh(T) 1 
sinh(T) n 2 
n 


a T E„- = - n 2 TT^ (^ + - —)s, 


sinh(T) 

ab 


(2NY,ij — Cxgij) 


sinh(T) 


N{Rij+ng ab - 2E ifc E^) (5.4) 


+ sinh(T) ^ Cx ^ ij n 9ij 


5.1.2 Rescaling of the reversed CMC flow 

In the reversed CMC case, we rescale with s(r) = 1 — (-) 2 . Then, some signs change. The constraints 
are 

R(g) — |E | 2 = (n — l)n, V a E ab = 0. (5.5) 

The defining equations for lapse and shift are 


A s iV=l + iV(|E| 2 -n), 

A g X i + K i m X m =2VjiVE i - 7 - sinh(T)(2 - n)V‘JV - ( 2NE mn 


(£ x g) mn )( r* 


r* 


(5.6) 
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and the evolution equations are 


9t 9ij = 
dr Sij = 


sinhCT) 
cosh(T) v 
2 sinh(T) 1 


cosh(T) 


(2 NUij — C x gij ), 


2 TV n 

cosh(T) n 2 n ab cosh(T) 

77 / 1 

eosMT) (£a ’ S+ “ V<Vj ' JV)- 


JV(i2y + n<7„ 6 - 2£ ifc E*) (5.7) 


5.2 Local existence 

We have the following local existence theorem in CMCSH gauge, for the initial data, which we consider 
in this paper. We distinguish between the cases of positive and negative curvature beginning with the 
latter. 

Lemma 5.1 (Analogous to I AnMo03 . Theorem 5.1]). Let 7 be a fixed Einstein metric on M such that 7 
is a metric of negative scalar curvature and s > n/2 + 1. Furthermore, let (go, kg) be CMCSH initial 
data on M such that 

|| flb - 7 || H . + ||E|| ff .- I <£ (5.8) 

with e sufficiently small. Then the CMCSH Cauchy problem is strongly locally well-posed in C k (FL s ), 
k = |_sj and the corresponding Lorentz metric g is a vacuum solution of the Einstein equations. The 
following continuation principle holds. There exists a S > 0 such that for [To,T+) being the maximal 
future existence inten’al to the given initial data at Tq in the rescaled time T, then either [To, 7+) = 
[arcsinh(l), 00 ) or 

limsup (||<7 - 7 ||jj/. + ||E||if»-i) >£ + S (5.9) 

for T —y T+ orT —>■ T_. 

The positive case is a bit more subtle due fact that the elliptic operators for lapse and shift and not 
necessarily isomorphisms. 

Lemma 5.2. Let 7 be a fixed Einstein metric on M such that 7 is a metric of positive scalar curvature, 
—2(n — 1) (f Spec(A 7 ) and 7 admits no Killing vector fields and let s > n/2 + 1. Furthermore, let 
(go, kg) be CMCSH initial data on M such that 

\\g — 7II-H'” + ll^l|i? 5 -i < £ (5.10) 

with e sufficiently small to assure that the conditions on 7 hold for all g € S|( 7 , 0). Then the CMCSH 
Cauchy problem is strongly locally well-posed in C k ('H s ), k = |_sj and the corresponding Lorentz metric 
g is a vacuum solution of the Einstein equations. The following continuation principle holds. There exists 
a S > 0 such that for (T_,T+) being the maximal existence interval to the given initial data at To in 
rescaled time T, either (T_, T+) = (— 00 , 00 ) or 

limsup (||<7 — 7 ||jj» + ||E|| ffs -i) > e + (5 (5.11) 


for T -> T + orT T_. 

Proof. The lemma for the negative case follows straightforward by the same methods as in lAnMo03ll . 
which would even yield a more general result without the smallness assumptions. 

In the case of positive curvature in the second lemma one has to assure that the elliptic operators defining 
the lapse and shift equation are in fact isomorphisms to use the relevant structure of the elliptic system as 
in lemma 3.2 of I AnMo03 1. The conditions we impose on the Einstein metric assure that these operators 
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are isomorphisms, c.f. section 2.5 As the perturbations are chosen to be small we can assure that the 


isomorphism property holds also for the perturbed data as long as we remain in an e-ball. This justifies 
the continuation criterion, which automatically covers the case implied by the analysis in HAnMo()3l . □ 


Remark 5.3. The local existence results mentioned above hold under more general conditions on the 
initial data (cf. I AnMo()3 l). We have stated a concise version which covers the case which is needed in 
the present paper. It is also understood that we choose the bootstrap assumptions in the proof of global 
existence such that the solution is contained inside the corresponding £-ball of the corresponding local 
existence criterion. 


5.3 Elliptic system 

We derive the elliptic estimates for lapse and shift in the case of positive and negative Ricci-curvature of 
the spatial metric. 

Lemma 5.4. Let s > n/2 + 2and (g, £) € B a s ( 7 ) x B s s ~ x (0) for some sufficiently small S g ,5s > 0, 
then 


ll-W — ~ll« < C(5 g , <Ss)||E||g_ 2 , 


1*11 


< C(Sg, Ss) 


12 VjNEV 


2NE mn (r mn - r m 


(5.12) 


■ ncosh(T)(l-)||V*iV|| s _ 2 


In combination with the first inequality this implies for X, 


l*|| a <<7(5 ff ,a E ) 


|£|IL 2 + 


,_2||E|| a _2||5-7ll-i 


+ ncosh(T)(l-)||£| 


(5.13) 


n 

In the case of positive curvature, cosh(T) is replaced by sinh(T) in the estimate. 

Proof We prove the estimates in the case of negative Ricci curvature. The positive case is analogous. 
The lapse equation, in the negative case, reads 

A 9 -/V = -1 + W(|E|g + n). (5.14) 

The maximum principle immediately yields the second estimate. Rewriting the lapse equation, we obtain 

A g (N -- n(N --) = N\E\ 2 g , (5.15) 

n n M 

which in combination with the point wise estimate on N and elliptic regularity for the operator A s — n 
yields the first estimate. 


Finally, we consider the estimate for the shift vector. We write the equation for the shift in the form 


A g A- + = F x + (C xg r n (r mn - r mn ), 

where 

Fx = 2X r \>y - ncosh(T)(l - 2 )V*iV - 2NT, mn (T i - f* m J. 

n 


(5.16) 

(5.17) 
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Elliptic regularity applied to the equation for the shift then implies 


< C{ 6 a ) [\\F x \\ H s-2 + \\(C x g) mn (F mn - V l mn )\\ H s-2 
Using the smallness of g — 7 we estimate 

\\{C x g) mn {F mn - f^JH ^-2 < C(6 g )\\X\\ H .-4g - 


(5.18) 

(5.19) 


By choosing the S g sufficiently small we can estimate the RHS by |||X||jj»-i and absorb it into the LHS 
above yielding 

\\X\\ H s<2C(S g )\\F x \\ Ha -2. (5.20) 

This finishes the proof. □ 


5.4 Energy estimate 

We restrict in the remainder to the case of negative curvature, the positive case is analogous. Before 
defining the total energy of the system, we cast the evolution equations into a form where the terms are 
ordered according to their eventual asymptotic behavior. This reads as follows. 


dr gij = 
dr = 


n cosh(T) , 

-2.ZVE,,- H-U-fA - 

sinh(T) J sinh(T) 


cosh(T) 

sinh(T) 

1 

sinh(T) 


(n — 1 )£ij + 


(CxFiij) + — 


sinh(T) 
n 


B, 


N (-n A g,j(g - 7 ) - R~i{g - 7 )) 


sinh(T) 
1 
2 


(5.21) 


c+ co ! MT) d _ 


sinh(T) sinh(T) 


Lemma 5.5. Let s > ^ + 1 and (g, £) £ 6 f ( 7 ) x BJ 1 (0 )for some S g ,Sy: >0 sufficiently small. Then 
for the perturbation terms, the following estimates hold. 


||A|| a <c(5 fl ,<y E ) 

»<2||VX|| S 

lie'll. 

PII.<e(<y fl ,<y E ) 


l|S||2- 2 (||ff-7ll. + ll7ll.) 

m 2 . + IWW|S||2 + ||S||2(||s - 7lU + Il7ll«) + \Ms\\9 7112+1 

'||S||2- 2 ||S||.' 


(5.22) 


Proof. The perturbation terms have the explicit form 


A = -2(1 - nN){g^ - 7 ^ + 7 ^) 

B = C x gij = \7 a Xb + V h Xa 

C = -ViVjN - 2NT, ik T,j + (N — - 7ij + 7ij ) + NJ ig 

D = -{n- 2){Nn- !)£„. 


Note that J is the perturbation term given in lemma 2.2 with the corresponding estimate ( |2.12| ). Now, the 
lemma follows from the estimates from lemma 15.41 . □ 


We define the main total energy. 
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Definition 5.6. Let s > n/ 2 + 1, then we denote 


s— 1 


k =0 


E a (fl, S) = || 5 - 7 ||! 2(7) + ^(-l) fe (S, A* 7 S) L 2 (5i7) + - £(-l ) k {g - 7, A* 7 (g - 7 ) W, t) 

= i + i + m. 


fc=i 


(5.24) 


Remark 5.7. Note that the energy is equivalent to the H s x H s 1 -norm of (g, E). 


Lemma 5.8. Let s > n/ 2 + 1 and ( g , E) £ H s x H s 1 be a solution to ( |5.3t > ( [5~4] >. Then there exists 
an e > 0 such that for 

(5,E)eBf( 7 )xSr 1 (0) J (5.25) 


the estimate 




(5.26) 


holds in the case of negative curvature 0/7. The analogous estimate with sinh(T) replaced by cosh(T) 
holds in the case of positive curvature of'y. 

Proof We consider the negative curvature case for the proof. The case of positive curvature is analogous. 
We take the time derivatives of the three individual terms. 

< 9 t ||9 - 7|||2 (7 ) = 2 [ (drg, (5 - 7)) 7 M7 

J M 

r n cosh(T) n „ . 

L + Sh(T )' A + is ~ 

j N(S, (9 - 7 )),99 + J (X (9 - 1))-,^ 


= 2 


m sinh(T) 
4 n 


(5.27) 


sinh(T) J M ' 

2n [ (B,(g ~ 

JM 


sinh(T) J M 

where A and B are given in (| 5.23 [>. In particular, we have an estimate of the form 


\dr\\9 - T\\h h ) \ < C 


1 


l~E s ( 5 ,E) + ^§||A|| i 2 ^E^E) 


Lsinh(T) ' sinh(T) 

1 '\B\\ L2 s/E s (g,X) 


sinh(T) 


(5.28) 


We order the terms according to their appearance in the energy estimate, substitute the expressions for A 
and B. 


I^tIIs - 7lll 2 ( 7) l < c 


1 


-nE s (y,E) 


sinh(T) 

cosh(T) ■ 


+ c '('Ms)—I|£|li 2 (lls-7|| S + Il7lls) v/E Jg^j 


9 ' sinh(T) L 
2 


sinh(T) 1 ~ v 9 ’ 


C(6 g ,6v) ||S||s _ 2 + ||iV|| s _ 2 ||E || s _ 2 || 5 - 7 |U-r 


(5.29) 


+ cosh(T)(l - - l)||E||g _ 3 ) ^E^E) 
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We sort the terms on the RHS into two different categories, where C is a new constant and we make use 
of the fact that E s (<?, £) < 1. Then we obtain 


\(h\\9 - 7lli,2( 7 )| < 


C 

sinh(T) 


nE s {g,X) + C 


cosh(T) 

sinh(T) 


l|S||WE s ( fl ,£) 


c cosh(T) 

sinh(T) 


||E||2_ 3 V / E^E). 


(5.30) 


Remark 5.9. The two terms on the RHS contribute either to the estimate to be proven or will be absorbed 
by a large negative term as demonstrated further below. 

We proceed with the evaluation of the next time derivative. 

s-l 


a T ^(-l) fc (S,A g fc 7 E ) L 2 (fl>7) 

k =0 




o / v' 


(5.31) 


f (N(A gn (g — 7)), Age'll) fj, g 
J M 


sinh(T) J M 

~ sinh(T) ~ 7)) ’ A s,7 E )^s + sinh ( T ) S M ^ Cx ^' A 9^^9 

+ sirdi(T) Jj C ’^ Jj °’ 

+ [ <£,[3r,A* 7 ]£)/i 9 +/ (£,A‘ 7 £)drM4 

7m 7m J 

Before we further evaluate the previous term, we compute the time derivative of the last term of the en¬ 
ergy and evaluate both terms in combination. 


The time derivative of the last term reads 


dT \“ 7 , A k gn {g - 7))l 2 ( 9 ,7) 
k=l 

= / 2 (^t( 5-7),A^ ( 5 i-7)) + ( fi (-7, [a T ,A* 7 ]( 5 -7))/r s 

2 fc =t l 

+ / ((5 — 7)) Ag j7 (<7 — 7))i9 T /i ff 1 
J M J 

- 5iRf) 

+ sin ^ (T) J m ( b ,Ag )7 (5f - 7 ))m 9 + J M (g - 7, [dr, Ag, 7 ](5 - 7 )}m 9 

+ [ ((5-7),A* 7 (p-7))ar/i ff l- 

7 m J 


(5.32) 
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In both previous computations there are commutator terms arising. We do an intermediate discussion of 
those terms in the following. The commutator operator can also be written as 

fc— 1 

Pr, A* 7 ](/i) = (d T A£ 7 )(/0 = X!(A' l7 o (5 t A s , 7 ) o A' 7 *- 1 )^) (5.33) 

1=0 

for some 2-tensor h. Recall that the Laplacian appearing here has the local formula 

A g^hij = —V[ 7 ]m(<7 m > s V[ 7 ]nM (5.34) 

L l g 

This shows that the variation of the operator with respect to the metric can be written schematically as 

(<9 t A S;7 )(/i) =d T g * A gn h 4- \/[^} m {g g d T g * V[ 7 ] n h) 

A 4 g (5.35) 

=&rg * V 2 [ 7 ]/i + d^g * V[ 7 ]g * V[ 7 ]/i + V[ 7 ]9t g * V[ 7 ]ft. 


where * is Hamilton’s notation of a combination of tensor products with contractions with respect to g. 
Therefore, by (|5.33|l and (| 5.35 [> and by suitable integration by parts, one can see that 


s — 1 

IE 

k=0 


(-1)" / <£,[3r,A* ]£> /ifl |<C'||£ 


IM 


I H* 


l (g) 




2 (g) 


(5.36) 


and similarly. 


iE( 

k=1 


(-!) fc / (5-7, [^t, A^ 7 ]( 5 - 7 ))/x 9 | < C 115-711^.^) \\drg\\ H s-i, g) (5.37) 

J M 


and for 8t g, we have good estimates which make this terms to be of higher order. 


Before we continue, we note the following estimate for the norm of the time derivative of g, which follows 
straightforward from the previous estimates. We have 




1 

sinh(T) 


|A'||ff S -i( g )||£||_y»-i( g ) 


c °sh (T) || E ||^ 


sinll ( T )..“ 1 .«*- i ( fl )( 1 + || s -' y || ^(»)) 

+ s iTlh (T) (ll S llg— Hfl) + ll JV 'llff- 1 (fl)ll s llff- 1 (fl)llff “ Tllff— 1(a)) 


(5.38) 


We have evaluated both commutator terms arising in ( |5.31| ) and ( |5.32[ > and proceed by combining the 
terms on the corresponding RHS into two different classes. We rearrange the terms. 
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dr(I + m) = — 2 


cosh(T) 

sinh(T) 


(n- l)X>l) fc / <E,A* 7 E)/x 

fc=0 *' M 


■ ^ (_1 )‘SS(T) /„ W(E ’ a ” (9 - 7) >'‘» 


s— 1 


+ E<-d‘ - 


k—0 


^mt)L { n{ 2^,(9 

- [ ((£,E),4*,E),, 

J M 


sinh(T) 7 m 
2 


/ (C, Ag i7 E)/ig 
JM 


sinh(T) 7 m 


+ [ (E,[a T ,A g fc i 7 ]E ) fig + f (E,A g fe 7 E)d TMg 

J M J M 


+ 5 D - 1 ) 1 


fc=l 


^cosh(T) 

J sinh(T) J M 

2 

+ 


[ (A A* ( 5 - 7 ))^ 

J M 

[ ( B , A gi 7 (5-7))M S 

JM 


< -2 


(5.39) 


sinh(T) 

+ [ (fl-7)[^r,A* 7 ](5-7))At g 

JM 

+ [ {(5-7)>Ag )7 (5-7))5r/i fl l 
J M J 

c °sh (T) (ri ^i)g(_i) fc f (E)A fe 7 E )Mg + C -^/E^E)||E ||| r ._ 1 

7 , n J M S --- / 


sinh(T) 


k—0 


(**) 


(*) 


c 


:E s (g, E) 


sinh(T) 

We still have to justify the last inequality. Therefore we have to analyze all terms on the right hand side 
above and estimate them by one of the three terms given in the last two lines of ( |5.39| >. Before we do this, 
we first argue, why this estimate implies the result. The first term on the right hand side, (*), has a nega¬ 
tive sign and its absolut value bounds the fT s-1 -norm of E up to a multiplicative positive constant from 
above. Therefore choosing i v /E s (g, E) sufficiently small by choosing the g close to 7 we can ensure that 
the second term, (**), is always bounded from above by the absolut value of the first term and thereby 
the sum of both terms is negative and can be estimated from above by 0 yielding the desired estimate. 

To complete the proof we need to justify the last estimate in (|5.39|>. We proceed term by term. 


The second and third line of the RHS of (|5.39|) contain leading terms of too high regularity to close the 
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estimate - these terms cancel pairwise using integration by parts. The resulting term can then be estimated 
as follows. 


sin h (T) f(l A 9,^9 - 7), N]A g ^(g - 7 ), E)/z s 


< 


C 


(5.40) 


sinh(T) “ TllffHfllllEllff-Hg) 

This term contributes to the last term on the RHS of ( |5.39| l. 

O 

The term in the fourth line is evaluated using R 7 (g — 7 ) 7 - = Rikjii^ig — "f) kl - We estimate 


-4 


sinh(T) 


S— 1 /* 

E (NR^(g- 7),a^ 7 s)/x £ 

fc— 1 J 


< 


c 


sinh(T) 


ff a (s)ll3 - 7l|rr s (a)||^||i? s - 1 (g)) (5.41) 


which contributes again to the last term on the right hand side of (|5.39|>. 


To evaluate the term in the fourth line one needs to observe the symmetry when using the integration by 
parts. This yields an estimate of the form 


2 

sinh(T) 


S — 1 p 

E / (M, A g fc 7 S ) M9 

k= 


< 


C 

sinh(T) 


ll*llff« (fl) E 8 (s,E). 


(5.42) 


The sixth and seventh line can be evaluated straightforwardly using the estimates for the H s norms of C 
and D as defined in ( |5.23| ). As the second term containing D does not contain a good time factor it is 
important to note that this term contains a factor quadratic in the fT s_1 -norm of E. Precisely, we have 


s-l 

^ sinh(T) J M 
2 




< 


cosh(T) 


s inh(T)^ C, ^ s " 1 (ff) H E ^ 5_1(s) + 2 sinh (T) H- P ll g<, - 1 (fl)H S H g,, - 1 (g) 


s slm l|E| 


(5.43) 




E s(g, S)(1 + ||iV||^ (ff) + y/E s (<?,£)) 


I o r COsh ( T ) ||v||4 

+ 2 C ^h(E 11 llff3 " 1(9) ’ 


which contributes to the second term and the third term on the RHS. Note, that by smallness of E we can 
absorb terms with higher exponents into the explicitly given one. 


The first term in the sixth line has been evaluated in ( |5.36| ) and clearly is quadratic in the /7 S_1 -norm of 
S. The other factors yield the necessary energy factor. 

We recall the estimate 


E 

fc =1 

which in combination with the estimate for the last factor, ( |5.38| > shows that the terms can be estimated 
as claimed contributing to the second and third term on the RHS of (|5.39|>. 


(E, [c*r, A 9 ]E)/r £ 


< C IIE 


Iff* 


l (g) 


11^11 


ff* 


2 (ff) 


(5.44) 
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The second term in the eigth line is determined by the time derivative of the volume form. We use the 
identity 


^="(Srti (JV + ^m v ‘ x 'h 


(5.45) 


and estimate 

|E / (s,A 9 fc 7 E)a TM£ 


k=1 


! M 




|E||^_ 1(ff) . (5.46) 


These terms can be estimated by the second and third term in (|5.39|l. 


The terms in the lines nine to twelve of ( |5.39| l are of the same type as the terms above. We briefly describe 
the way to estimate them in the following. 

The ninth and tenth line can be estimated as follows 


J sinh(T) 

< C 


-sinh(T) 


Ij A - O - + SS(T) Ij B • A ‘> - 

B \\H°(g)\\g - 7l|ff s (g) 


COSh(T) P||^ (g) || 5 - 7 ||^ (g) + ^ 1 


(5.47) 


sinh(T) 


In combination with the estimates for A and B we can estimate these terms by the terms on the right 
hand side of ( |5.39| >. The commutator term in the eleventh line has been already evaluated in ( |5.37[ ), which 
yields an estimate of the form as treated above. Finally, the last line simply contains the time derivative 
of the metric and can be treated as the corresponding term above. 


We have analyzed all relevant terms in the estimate, which in combination with the argument following 
(|5.39|> finishes the proof. 


□ 


5.5 Improved decay 


We proceed by deriving an energy estimate for the Sobolev norm of the trace free part of the second 
fundamental form in one order of regularity below the maximal regularity. This estimate holds under 
the condition that boundedness of the total energy of maximal regularity is given. The structure of this 
estimate is such that it eventually leads to decay of the trace free part of the second fundamental form in 
this regularity. 


Lemma 5.10. Let 


Assume for some 0 < £ < 1 


5—2 


H s - 2 (E) = ^(-l) fe (E, A£ 7 E) L 2 (fli7) . 

k=0 

E s (5,E) < £. 


Then 


< —2(n - l)ggH._ 2 (E) + egg 
C 


(h s _ 2 (E) 


-^^||^i (s)Hs _ 2 (E) 

+ (1 + ||-^'||fl’«-l(g) + v/^)v / ^'\/Hs- 2 (S) 


(5.48) 

(5.49) 


(5.50) 
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Proof. The proof follows from the identical computations in the proof of lemma 5.8 The higher order 
term containing the Laplacian is estimated using integration by parts and Holder’s estimate. □ 

In combination with the elliptic estimates for lapse and shift ( |5.4| > we obtain the following corollary. 

Corollary 5.11. Under the assumptions of lemma\5. 10 the following estimate holds. 


9 t H s _ 2 (E) < —2(n — 1) 


cosh(T) 

sinh(T) 


H S _ 2 (E)+C 


cosh(T) / 
sinh(T) V 


H 


s -2 


(£))' 


+ 


c 


sinh(T) . 


(H S _ 2 (E) + yfe(C + H s _ 2 (E))) (h s _ 2 (E ))' V2 (5-51) 


As in particular 


+ (1 + H s _ 3 (E) + v / e)Ve\/H s - 2 (E) 
H s _ 2 (E) < e 


this simplifies to 

In the following we deduce an estimate for H S _ 2 (E). 

Lemma 5.12. Under the same assumptions as in the previous lemmas and if n > 2, 


(5.52) 


(5.53) 


\/H s _ 2 (E) < 


If n = 2, the estimate is 


v/H s _ 2 (E) < 


C y/s 
sinh(T) ‘ 


Cy/e 


sin. 


h 1 / 2 (T)' 


Proof. If n > 2, we have 


9r AWSj <-S^V^^(5j + 


and if n = 2 , we have 




(5.54) 


(5.55) 


(5.56) 


(5.57) 


provided that e is small enough. These differential inequalities immediately imply the desired decay. □ 
Lemma 5.13. Under the same assumptions as above, the shift vector admits the estimates 

||X|| g < Ce ■ sinh _ 1 (T) (5.58) 

if n> 2 and 

\\X\\ e < Ce-sinh- 1 / 2 ^), (5.59) 

ifn = 2 . 

Proof. This follows from the previous lemma and ( |5. 1 3[ >. □ 

Remark 5.14. From the decay of X and E, we also get an exponential decay of OtH in the // s ~ 2 -norm, 
i.e. the rescaled metrics converge to a limit metric in H s ~ 2 as T — > oo. This follows from ( |5.21| > and 

1(532}. 
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5.6 Proof of the main theorem 


Using the previous lemmas we are now able to state the proof of Theorem 1 1.1 1 and Theorem 1.2 We give 
the explicit proof for Theorem |1.1| the positive case is analogous. 


Before we consider r)-small CMCSH-initial data at initial 

s'-l 


2.9 


we have a small H s x H s 


Proof of theorem \L1\ Let e > 0 be fixed. 

time v/2, assume we start with arbitrary initial data. By Theorem 
neighbourhood V in the set of arbitrary initial data such that the MGHD of any data in V admits a 
hypersurface of constant mean curvature — sfl and the induced data {g 0 , E 0 ) stays in a small H s x H s ~ l 
neighbourhood U of the initial data of the background. By Theorem |2. 5 1 we can pull back the data along 
a diffeomorphism p such that (p*go, p*T,o) satisfies the CMCSH gauge and is <5-close to the initial 
data of the background solution. From now on, we let the data evolve under the Einstein-flow and the 
corresponding solution will be isometric to the MGHD of the initial data we started with. 

Without loss of generality, assume that e is so small that ( ]5.26| ) holds as long as E a (g(t), S(f)) < 2e. 
Let T max be the maximal existence time of the solution and suppose that T max < oo. By the Gronwall 
inequality, 

£ s (<?(T max ),£(T max )) < e c(2e)/ ^ max ™> . Es (g(V 2),E(V2)) 

(j.oU) 

<e c(2 e)f£ sEEfry . E s (g(V 2), £(V2)) 

which shows that the left hand side is bounded by some arbitrarily small e\ supposed that S was chosen 
small enough. Due to local existence, this contradicts the maximality of the existence time. Therefore, 
T max = oo and since the energy is equivalent to the H s x // s_ 1 - norm of the data, we also obtained the 
desired bound on the solution. 

To complete the proof we show global hyperbolicity and future- and null geodesic completeness. For 
this purpose, some properties of the lapse, the shift, the second fundamental form and the family of 
Riemannian metrics have to be checked. At first, by the estimate of the main energy, lemma |5.4| and 
lemma |5T3| we can choose for any e > 0 a neighbourhood of the initial data such that 


Il5-7llc°( 7 ) < £ : 


N — 


1 


< £sinh *(T), ||ViV|| c o( s) < esinh 1 (T), 


n c » 

llc»(g) ^ £ sinh” 1 / 2 (T), ||X|| C0 ( S) < £sinfr 1 / 2 (T) 


(5.61) 


for all T > T 0 . Recall that the non-rescaled metrics g satisfy g = s(r) 1 g = sinh 2 (T)g which implies 
that for all tangent vectors v, we have 


7 ijV l v : * < CgijV l vi < C sinh 2 (T )g i jV z v :i 


(5.62) 


which implies that the metrics g are uniformly bounded from below by 7 . For E = s 1 / 2 (r)E, we have 

|E| S < sinh _ 1 (T) • |E| S < sinh" 3 / 2 (T) • £ (5.63) 


For the original shift vector X = sinh(T)X, we have 


X 


C°(g) 


< esinh 3 / 2 (T). 


(5.64) 


However, this is the shift vector from the CMC-gauge and we have to define the shift vector which 
corresponds to the time function T: 


X 9 :=( -L,.) = 


/ d 


<9T’7 sinh 2 (T) \5t’ 7 sinh 2 (T) 


-XL 


(5.65) 
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Here, and X' 1 are the 1-forms which are equivalent to X and X via the metric g. Now we immediately 
get 


||X||. < en ■ sinh 1 ^ 2 (T) 

For the original Lapse function N we have N = s{t)~ 1 N = sinh 2 (T)iV, which yields 


sinh 2 (T)7V — — 

< esinh 1 (T), 

VN 

n 

c° 



C°(g) 


< £. 


(5.66) 


(5.67) 


However, as for the shift vector, N was obtained by the CMC- gauge. We now compute the lapse-function 
N according to the time T. We have 




n 

xx. 


,sinh 2 (T) 

which immediately implies 

||N — 11 |< e ■ nsinh _1 (T), 




(5.68) 


V7V 


<e-nsinh 2 (T). 


c°(g) 


(5.69) 


By [ ChCo02 , Theorem 2.1], global hyperbolicity follows from ( |5.62 1 , ( 5.69| ) and ( | 5.66 [ ). Causal com¬ 
pleteness of the solutions follows from ( |5.63| > and ( |5 ,69[ > by using !!ChCo021 Theorem 3.2 and Corollary 
3.3], □ 
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